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ABSTRACT 

Type la supernovae (SNe la) play a crucial role as standardizable cosmological candles, though the 
nature of their progenitors is a subject of active investigation. Recent observational and theoretical 
work has pointed to merging white dwarf binaries, referred to as the double-degenerate channel, as 
the possible progenitor systems for some SNe la. Additionally, recent theoretical work suggests that 
mergers which fail to detonate may produce magnetized, rapidly-rotating white dwarfs. In this paper, 
we present the first multidimensional simulations of the post-merger evolution of white dwarf binaries 
to include the effect of the magnetic field. In these systems, the two white dwarfs complete a final 
merger on a dynamical timescale, and are tidally disrupted, producing a rapidly-rotating white dwarf 
merger surrounded by a hot corona and a thick, differentially-rotating disk. The disk is strongly 
susceptible to the magnetorotational instability (MRI), and we demonstrate that this leads to the 
rapid growth of an initially dynamically weak magnetic field in the disk, the spin-down of the white 
dwarf merger, and to the subsequent central ignition of the white dwarf merger. Additionally, these 
magnetized models exhibit new features not present in prior hydrodynamic studies of white dwarf 
mergers, including the development of MRI turbulence in the hot disk, magnetized outflows carrying 
a significant fraction of the disk mass, and the magnetization of the white dwarf merger to field 
strengths ~ 2 x 10* G. We discuss the impact of our findings on the origin and observed properties of 
SNe la and magnetized white dwarfs. 

Subject headings: supernovae: general — supernovae: individual (2011fe) — magnetohydrodynamics 
— white dwarfs — ISM: supernova remnants 



1. INTRODUCTION 

Type la supernovae (SNe la) are among the most ener- 
getic explosions in the known universe, releasing 10^"'^ erg 
of kinetic energy and synthesizing 0.7 of radioactive 
^^Ni in the ejecta of a typical bri ghtness explos ion. The 
discovery of the Phillips relation ( Phillips||1993 1 enabled 
the use of Type la supernova (SIN la) as standardizable 
cosmological candles, and has ushered in a new era of 
astronomy le ading to the discovery of the acce leration of 
the universe ( |Riess et al]|1998[ [Perlmutter et al 1999). 

The nature ot the 'i'ype la progenitors, as well as their 
precise explosion mechanism, remains a subject of active 
investigation, both observationally as well as theoreti- 
cally. Observational progress to determine the nature 
of the SNe la progenitors, as well as their underlying 
explosion mechanism, has accelerated in recent years, 
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with a series of projects, including the Palomar Transient 
Factory (PTE), Large Synoptic Survey Telescope, Pan- 
STARRS, and the Dark Energy Survey, all coming on- 
line. This progress culminated in th e discovery of 2011f e 
in MlOl by PTF on August 24, 2011 ( [Nugent et al.|2011[ ). 
At a distance of 6.4 Mpc, 2011fe is the nearest SNe la 
detected in the last 25 years, and has proven to be the 
kind of SNe la exemplar system that SN 1987A has been 
for SNe II. PTF captured 2011fe within 11 hours of the 
explosion, making it the earliest SN la ever detected, 
and opening the gates to prompt multi-waveband follow- 
up observations in the radio, optical, UV, and X-ray 
bands. The first weeks of multi-wavelength follow-up ob- 
servations have directly confirmed for the first time that 
the primary object is a carbon-oxygen white dwarf, and 
have placed tight constraints on the progenitor system to 
201 Ife, ruling out red giant as well as Roche- lobe over 
flowing main seque nce companion s (Nug ent et al 
BrownlTt al. 2012^ jHoresh et al. 2012 ,Li et al. 
Bloom et al.|2012| )' 

Additionally, determinations of the SNe la delay time 
distribution (DTD) generally follow a. power-law, 
consistent with expectations for double-d egenerate (DD) 
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systems fc al-Yam & Maoz 2004; Tot ani et al.||2008 



Maoz & B adenes||2010| |Maoz et al.||201o| ). A relate? 



but independent model ol the supernovae rate based 
upon a two-component model accounting for both a 
prompt and delayed component a lso supports the exis 
tence of a delayed, DD channe l (Scannapieco & Bild 

Moreover 



sten 
lor 



2005[ Raskin et al.||2009a [). Moreover, the search 
30th point sources m pre-explosion archival data 
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and "ex-companions" in SNe la remnants have so far 
yielded no definitive candidates, wit h the possible con- 
tested example of Tycho's r emnant fMaoz fc Ma nnucci 
2008; Gonzale z Hernandez et al. 2012; Schacfcr fc .FaK- 



notta 2012t |Edwards et al.|2012(|Kerz endorf ct al.|2012p . 



Tastly, the search tor hydrogen m the nebular spectra of 
remnants places fairly tight constraints on the amount of 
hydroge n (< 10"^ Mp ) than can be stripped from a com 



panion ( Leonard 2007 1 . It may be possible to understand 
both the absence ot ex-companions and nebular hydrogen 
in the context of the single-degenerate channel as due to 
the time delay of the spin-down of th e white dwarf after 



accretion from its companion ceases (Di Stefano & Kilic 



2012 1. However, the weight of the observational evidence 



strongly suggests the viability of the DD channel model 
as the progenitors for some, if not the majority, of SNe la 
events. 

The key conceptual challenge faced by the DD channel 
for SNe la is to explain how these models yield a ther- 
monuclear runaway, as opposed to an accretion-induced 
collapse (AIC) to a neutron star. Early spherically- 
symmetric models, based on Eddington accretion rates 
onto the white dwarf merger, suggested that double de- 
generate mergers will ignite a carbon-burning deflagra- 
tion wave which propagates inward to the core of a 
1.0 Mq 50/50 carbon-oxygen white dwarf in ^ 2 x 10* yr, 
comparable to the thermal tim escale of the white dwarf 
merg er, resulting in a AIC 



( S aio & Nomotol 
An AIC may be avoided, however 
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2004f. 

evolution within the rotating merger and fully multi 
dimensional, magnetized accretion disk differs signifi- 
cantly from the one-dimensional models, and can ignite a 
detonation on a timescale much shorter than the inward 



carbon deflagration timescale (Nomoto & Iben 1985 1 . An 
additional crucial difference highlighted by more recent 
numerical simulations deals with the thermal structure 
of the white dwar f merger itself 
19891 [T990l|Ras io 
Guerrero et al 



2007; 




Shapiro|[T995 
20041 iD'Souza et al 



^007; 
2011). 



Mochkovi tch fc Livio 
'grctain et al.||1997 
120061 
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Pakmo r et a] 

These simulations demonstrate 



Yoon et al 
Dan et al. 



20111 |Zhu et a 
that the comprcssional work produced during the final 
tidal disruption of the white dwarf binary results in a 
hot (~ 10® K) rotating white dwarf merger, quite unlike 
the cold (lO'' K) isothermal white dwarfs typically taken 
as a starting point in earlier one-dimensional studies. 

Furthermore, recent work suggests that the outcome of 
white dwarf mergers may not always be either a SNe la 
or an AIC, but could also result in a high-field mag- 
netic white dwarf (HFMWDs). HFMWDs have magnetic 
fields in excess of 10^ G and up to 10^ G. Very few of these 
white dwarfs belong to a non-interacting binary system, 
and moreover th ey are more massive than average 



see, for instance, Kawka et al. (2007). All these char- 
acteristics point towards a binary o rigin of these white 
dwarfs. Although long-suspected ( Wickramasinghe & 



Ferrario|[20 00), it has only recently been shown that it 



the white dwarf merger fails either to detonate into an 
SNe la or collapse down into an AIC, it will result in a 



magnetized , rapidly-rotating white dwarf ( Garcia-Berro 
et al.||2012 ). Whether the magnetic white dwart rotates 
rapidly or not depends on the relative orientation of the 
magnetic and rotational well as on the efficiency 



of the various braking mechanisms. 

The evolution of the white dwarf merger subsequent 
to the coalescence of the initial binary system remains a 
subject of active investigation. Numerical models have 
begun to relax assumptions of earlier work by model- 
ing the accretion of the hot thick accretion disk onto the 
white dwarf, either by including a prescription for t he ac- 
cretion proc ess and the spin-down of the merger (Yoon 
et al. 2007), or by employing a Shakura-Sunyaev tur - 
bulent viscosity (Shen et al. 2012 Schwab et al. 2012). 
Other researchers have investigated the violent rne rgers 
of super-Chandrasek har mass white dwarf systems ( Pak - 
mor etal.|20lol|2011[ ) and colhsions ( Raskin et al.|200^Ff : 
While violent mergers and collisions ot white dwarts are 
found to lead to detonations, sub-Chandrasekhar merg- 
ers are more prevalent in nature. Recently, van Ker- 
jwijk and colleagues have re-invigorated the examina- 
tion of whether the accretion of the disk may give rise 
to a detonation, even for sub-Chandrasekhar mergers 



( van Kerkwijk et al.j|20To{ [Zhu et al.|2011 [van Kerkwijk| 
z012| ). Specifically, beginning with a near-equal mass bi- 
nary with two 0.6 Mq carbon-oxygen white dwarfs, which 
both subsequently tidally disrupt and merge, van Kerk- 
wijk et al suggest that the accretion of the thick, turbu- 
lent disk surrounding the white dwarf merger results in 
the comprcssional heating of the degenerate material in 
the white dwarf over a viscous timescale, which in turn 
leads to a detonation. The model explains how a realis- 
tic multidimensional DD merger might produce a SNe la 
instead of an AIC. Equally interesting is the possibility 
that sub-Chandrasekhar DD mergers may help to bring 
observed and predicted SNe rates in cl oser agreement 
( [Ruiter et al.|2009[|Ba"denes fc Maoz|2012[ ). Furthermore, 
detonations of sub-Ghandrasekhar progenitors naturally 
produce nucleosynthetic yields and luminosities closely 
in line with observation, and thereby sidestep the long- 
standing problem of pre-expansion required during the 
deflagration phase of near-Chandrasekhar mass white 
dwarf pr ogenitors encou ntered in the single-degenerate 
channel ([Sim et al.||2010|). 

Despite these advances, all previous numerical simula- 
tions to-date on double-degenerates have treated the in- 
fluence of magnetic field only approximately through the 
use of a Shukura-Sunyaev a prescription, or neglected its 
contribution entirely. Consequently, fundamental ques- 
tions involving the influence of the magnetic held upon 
the outcome of the merger remain unresolved. These 
questions include: What is the structure of the mag- 
netic field in the white dwarf merger, disk, and corona 
— ordered or disordered? If ordered magnetic fields are 
present, are they capable of coUimating outflows? Are 
significant amounts of mass outfluxed from the binary 
system? How does the magnetic field influence the spin 
of the merger? What is the influence of the field upon 
nuclear burning? Under what conditions may we gener- 
ally expect mergers to produce stable HFMWDs, SNe la, 
R Cor Bor stars, or accretion induced collapses? In this 
paper, we present the first set of numerical simulations 
of the post-merger evolution of double-degenerates to in- 
clude the effect of the magnetic field, and in so doing, to 
begin to address the role of the magnetic field upon each 
of these fundamental questions. 
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2. SIMULATIONS OF THE DOUBLE-DEGENERATE 
MODEL 

2.1. Simulation Setup 

The astrophysical fluid framework code FLASH has 
previously been used to simulate single-degenerate mod- 
els of SNe la in both 2D cylindrical a nd 3D Cartesian 
geometry in a wid e rang e of studies dTow nsley et ah 
20071 IJordan et al.lf2008; Mcakin ct al. 2009;T<'ai'ta et al. 



dB 

~dt 



V • {vB - Bv) = 



(4) 



Jordan et al.|2012a,bj . Here we build upon and ex- 
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tend this body of work to simulate the double-degenerate 
channel. Our double-degenerate models take as an ini- 
tial condition the endpoint of a 0.6 Mq + 0.6 Mq 40/60 
carbon-oxy gen white dwarf merge r from the 3D SPH sim- 
ulations of Loren-Aguilar et al.| ( 2009 ). The stars are 
modeled using 4x 10° particles. The initial condition cor- 
responds to non-sychronously rotating white dwarfs on 
a circular orbit with an orbital separation of ~ 0.03 i?©, 
which immediately leads to an unstable mass transfer 
episode lasting for ~ 500 s. The final configuration of the 
merger corresponds to a solidly-rotating central compact 
object (O ~ 0.3 s^^) surrounded by a thick accretion 
disk extending up to a distance of ~ 0.07 i?©. Unlike 
unequal-mass and synchronously- rotating merger events, 
the remnant of the irrotational 0.6 + 0.6 Mq case presents 
a temperature peak (~ 6 x 10^ K) at the center of the 
central compact object, with a rapid drop towards the 
outer parts of the disk. Whether or not the initial white 
dwarf binary is synchronously- rotating depends upon the 
efficiency of tidal dissipation, which remains an open 
question. Because the initial temperature profile plays 
a crucial role in determining the nuclear burning within 
the white dwarf merger, the possible influence of initial 
synchronous rotation of the white dwarf binary is an im- 
portant assumption of our models which must be borne 
in mind. No noticeable nuclear processing was found 
during the merger process. 

We utilize the SPH smoothing kernel to interpolate the 
Lagrangian SPH data onto a 2D axisymmetric cylindri- 
cal coordinate (r, z) Eulerian mesh, averaging over az- 
imuthal angle 0, while retaining all three velocity com- 
ponents {vr, W0, Vz) — see appendix for details. Both 
the interpolation and the azimuthal angle-averaging are 
not guaranteed to conserve energy; however, we have 
confirmed that the initial total energy on the mesh is 
preserved to within 1% or better of the SPH value for 
all models presented here. We then advance this ini- 
tial Eulerian initial condition in time using the adaptive 
mesh refinement (A MR) FLASH appl ication framework 
pubey et al.l200'9]|Fryxell et al.|2000| using an initially- 
weak poloidal magnetic held to seed the growth of the 
MRI. 

We solve the fundamental governing equations of self- 
gravitating, inviscid ideal magnetohydrodynamics, which 
can be expressed as: 



Here, p is mass density, v is the fluid velocity, B is the 
magnetic fleld, and g is the gravitational acceleration. 

= p + B"^ /{Sir) is the total pressure, including both 
gas pressure p and magnetic pressure / (Stt). pE is the 
total energy density, E = l/2pv^ + pe + B'^/{%tt). The 
inclusion of Poisson's equation for self-gravity: 



(5) 



where g — — V(/), as well as an equation of state close the 
system. Equation ([4]), the induction equation, preserves 
V • B = if this is imposed as an initial condition. 

The poloidal magnetic flcld is cfflcicntly initialized in a 
divergence-free form by dcflning the toroidal component 
of the vector potential A: 



Mr) 



BQ{p-po)f{r) ifp>po 
ifp</5o- 



(6) 



This form is motivated by previous MRI studies (Haw- 
ley 2000), with the inclusion of a fllter function f\r) = 
A tanh [(r — rg) /A]", chosen to localize the initial 
poloidal fleld to the disk, and to avoid initially magnetiz- 
ing the merger, the axial region near the z-axis, as well as 
low-density regions outside the disk. Here we have chosen 
ro — 10 km, A = 1.5 x 10^ km, po — 2 g/cm^, a = 9, and 
Bq is an overall fleld strength factor chosen to ensure the 
magnetic pressure is everywhere weak compared to the 
gas pressure initially. The magnetic fleld is then straight- 
forwardly defined at cell edges by finite-differencing the 
vector potential using B = V x A, and advanced using 
the unsplit id eal MHD solver in angular-momentum con- 
servmg form (Lee & Deanc 2009) [Tzeferacos et al ||2012] 
|Leej|2012j). The Roe Riemann solver is employed, with 
piecewise parabolic (PPM) spatial reconstruction and a 
minmod slope limiter. The divergence-free prolongation 
of the magnetic fleld is done using an a dapted implemen- 
tation of the method oflLi & Lil (|2004|). 



dpv 



dpE 



Our simulations employ an equation of state that in- 
cludes contributions from blackbody radiation, i ons, and 
electrons of an arbitrary degree of degeneracy (Timmes 
& Swesty 20001, along with an axisymmetric niultipole 
treatment ot gravity, with the series truncated after ten 
moments {(, = 10). Nuclear burning is incorporated 
through the use of a simplifled 13-species alpha-chain 
network, which includes the effect of neutrino cooling 
( |Timmes| p99). 

We performed a series of simulations from the 
previously-described initial condition, varying both our 
choice of the initially-weak magnetic fleld as well as 
the spatial resolution. The full set of completed runs 
is shown in table [l] All runs are performed on a do- 
main r < 1.31 X 10^*^ cm in the radial direction, and 
—6.55 X 10^ cm < z < 6.55 x 10^ cm in the vertical, with 
diode boundary conditions, which strictly guarantee that 
no inflow occurs at the outer boundary. The white dwarf 
merger has a radius ~ 1.5 x 10^ km, which is resolved 
with 60 cells in our standard model. The isothermal disk 
scale height H = ^/2cs/^ ~ 4.8 x 10^ cm at a temper- 
ature of 5 X 10^ K near the midpoint of the disk, where 
V • [v{pE+p^) - B{v ■ B)] ^ pg ■ V (3) ~ 0.05 s"^ At an inner disk radius of 2 x 10^ cm, the 



dp 
dt 



V • {pv) = 



V • {pvv - BB) + Vp* = pg 



(1) 



(2) 
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TABLE 1 
Simulations Properties 



Simulation Resolution ^ Ax(km)'' /J*^ i?mag (erg) H/Ax ° max (Ac/ Ax) ' Duration (s) ^ 



S-bh 


512^ 


256 


2000 


2.79 X lO"*"^ 


18.7 


76.3 


2 X 10* 


H-bh 


10242 


128 


2000 


2.79 X 10*6 


37.4 


158.2 


5 X 10^ 


L-bh 


2562 


512 


2000 


2.79 X 10*6 


9.4 


38.3 


5 X 10^ 


S-bm 


512^ 


256 


1000 


5.58 X 10*6 


18.7 


107.9 


10* 


S-bl 


512^ 


256 


500 


1.12 X 10*^ 


18.7 


152.6 


1.5 X 10* 



^ Number of grid cells per domain. 
^ Linear spatial resolution. 

^ Initial global average of the dimensionless ratio of gas pressure to magnetic pressure, over the entire spatial domain. 
^ Initial total magnetic energy on the problem domain. 

^ Ratio of isothermal disk scale height H to Aa;, evaluated near middle of disk. 
^ Initial maximum ratio of the effective resolution Ac/Aa:. 
^ Maximum duration of run. 



Keplerian period around a 1.0 Mq white dwarf merger is 
51 s. The viscous accretion time, based upon a simple 
constant Shakura-Sunyaev a = 0.01, is roughly 7 x 10'^ s. 
Our standard model (designated "S-bh") has a resolution 
of 256 km and an initial value of the ratio of gas pressure 
to magnetic pressure, defined to be the dimensionless 
ratio /3 = 8ttP/B^, of 2000. /3 varies throughout the spa- 
tial domain, but the magnetic pressure is initially signif- 
icantly less than the initial total gas pressure in all mod- 
els; its initial minimum value is 16.5 in model S-bh. We 
note that while our seed magnetic field is initially dynam- 
ically weak everywhere, and its initial magnitude in the 
white dwarf merger (2.8 x 10^ G) is typical of field white 
dwarfs (< 10^ G), its value in the disk is astrophysically 
large by this same comparison. Our choice is motivated 
by the requirement to accurately capture the dynamics of 
the MRI by well-resolving the fastest-growing MRI mode 
Ac — 6.49 va/^, where va is the local Alfven speed , 
and the local rotational velocity ( Hawley et al.|[l995 ), 
on a computationally-tractable mesh size. MKl simula- 
tions which do not initially resolve Ac also become unsta- 
ble and reach magnetic field saturation, but take much 
longer to do so, since only a narrow band of all unsta- 
ble modes, not includ ing the fastest-grow ing mode, are 
captured on the mesh (Suzuki & Inutsuka 2009). In real- 
ity, because the fastest-growing mode ot the MRI always 
grows on a dynamical timescale ^ we expect that 

even an astrophysically-realistic magnetic field strength 
will reach saturation on a relatively short timescale in 
comparison to the viscous timescale. 

Our additional models vary both the spatial resolution 
and the initial magnetic field strength, in order to test 
spatial convergence, as well as sensitivity to the resolu- 
tion of Ac- Our standard model S-bh maximally resolves 
Ac with approximately 76 cells per wavelength, and the 
disk scale height H (evaluated near the midpoint of the 
disk) with roughly 19 cells. The vertical resolution of 
our standard model is therefore somewhat less than 3D 
convergence studies of stratified shearing-box MRI simu- 
lations, which demonstrate that between 32 and 64 cells 



per sc ale height are required for convergence ( Davis et al 
2010 1. Moreover, because both the disk scale height in a 



global disk geometry and the initial seed magnetic field 
vary, we do not necessarily resolve either the scale height 
or the fastest-growing mode uniformly throughout the 
domain. To address this issue, we examine our conver- 
gence of the peak magnetic field and magnetic stresses by 



varying our resolution explicitly, increasing and decreas- 
ing our standard model resolutions in runs H-bh and L- 
bh, respectively. Additionally, we also vary our initial 
magnetic field strength, increasing it in both models S- 
bm and S-bl, which have initial /3 values of 1000 and 
500, respectively. Because Ac depends on the initial field 
strength, these model variations also vary the effective 
resolution (see table fll column f). 

Our standard model S-bh has been advanced to 2 x 
10** s, while other models run for varying durations. Our 
standard model duration is equivalent to 390 inner rota- 
tional periods, and several viscous accretion timescales. 
We output the state of the system in 10 s intervals. This 
extensive time series is sufficiently long to permit accu- 
rate time-averages over turbulent quantities. We exper- 
imented with both relaxed and "cold-start" initial con- 
ditions; due to the rapid growth of the MRI, there were 
relatively small changes in the outcome between these 
initialization procedures. The results presented here are 
all cold starts. 

2.2. Simulation Results 

2.2.1. Magnetic Structure of White Dwarf Merger and Disk 

A snapshot depicting the logarithm of the density 
taken from the midpoint of our standard model, at 
t = lO'' s, is shown in figure l(a)| whereas the distri- 



bution of temperatures is shown m figure |l(b) The 



rotating, hot white dwarf merger is surrounded by a 
differentially-rotati ng, t hick accr etion disk. On the bot- 
tom panels, figures 1(c) and 1(d) we reveal the magnetic 
structure of the merger and the disk. This is done plot- 
ting the poloidal magnetic field lines in the r-z plane 
superimposed on the background toroidal magnetic field 
Bcf, (left-hand panel), and the ratio of the gas to mag- 
netic pressure (3 (right-hand panel). Regions of high 
/3 are dominated by gas pressure, while those with low 
low P are supported by magnetic pressure. The merger 
and accretion disk themselves remain relatively weakly- 
magnetized (/3 ^ 1), whereas the disk corona and bi- 
conical jets are strongly- magnetized (0 < 1), a s previous 



MRI studies have found (Miller & Stone 



2000) 



The magnetic field structure in the disk is highly- 
turbulent and disordered. Loops of low-density, heated 
magnetic fl ux rise buoyantly abov e the accretion disk into 
the corona (Machida et al. 2000), where some reconnect 
through numerical resistivity, thereby heating the coro- 
nal region. Some poloidal loops of flux — which actually 



5 



le9 



le9 




0.6 

r (cm) 

(a) logp 



0.6 0.1 

r (cm) 

(b) logT 



le9 




le9 



0.0 



0.2 



0.4 



0.6 0.8 

r (cm) 

(c) Magnetic Field 



1.0 



1.2 



lelO 




0.6 0.8 

r (cm) 

(d) log/3 



Fig. 1. — Four frames in the r-z plane consisting of a) logp, b) logT, c) magnetic field, with lines of poloidal magnetic field in the r-z 
plane superposed against a color raster plot of the toroidal field B^, and d) the ratio of gas pressure to magnetic pressure /3 value. All four 
frames are taken at the midpoint of the model S-bh simulation at t = 10^ s. 



are toroidal in shape in an axisynimetric geometry — are 
long-lived in our simulation, persisting for many local dy- 
namical times. While it is known that these poloidal flux 
tori are subject to a wide variety of instabilities in 3D, 
including the kink and interchange instabilities, b oth the 



toroidal fi eld and the differential shear in the disk ( Spruit 
et al.||f995) may help stabilize these even in full SIT 



in contrast, biconical axial outflows carry open field 
lines away from the merger. The biconical region is 
strongly magnetized and heat ed to T ~ 10® K, as is 
clearly seen in figures 1(b) and 1(d) A strong outflow is 



driven at the interface of this region with the magnetized 
corona, similar to previous MRI stud ies of black hole ac- 
cretion disks ( jPe Villiers et al^|2005 ). However, we find 
that this interface region, which is Kelvin-Helmholtz un- 
stable, varies significantly in location and shape over the 
duration of the simulation. Moreover, although there is 
a net outflux of both mass and magnetic flux from the 
simulation domain, there are also thin returning flows, 
driven by reconnection, which infall onto the disk and 
merger, similar to that seen in previous work (Igumen- 
shchev et alT]|2003[ ). 
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Fig. 2. — A schematic diagram illustrating the key features of our magnetized model, after the development of the MRI. A weakly- 
magnetized Keplerian disk surrounds the uniformly-rotating white dwarf merger. The development of the MRI within the disk gives rise 
to strongly-magnetized biconical jets and a magnetized corona. Additionally, strong outflows are driven at the interface of the jet and the 
corona. 




Fig. 3. — a) The evolution of the global magnetic energy for all models (S-bh black, H-bh blue, L-bh green, S-bl cyan, S-bm red). The left 
panel inset shows the evolution of the magnetic energy for our standard model over the entire duration of the simulation, b) The effective 
Shakura-Sunyaev magnetic alpha coefficient am (see text for definition) for all runs, with the same color notation as the panel a. The right 
panel inset shows the evolution of am over the entire simulation for our standard model. 



Throughout this paper, we will have the need to 
separate the domain into the white dwarf merger, 
rotationally-supported disk, and magnetized corona in 
order to determine the properties of each of these re- 
gions individually — see figure 2. We characterize each 
of these regions based upon a physical criterion: specifi- 
cally, we define our magnetized coronal region to be dom- 
inated by magnetic pressure (/3 < 1), with the remainder 
of the weakly-magnetized gas {/3 > 1) divided into both 



the white dwarf merger and the disk. The white dwarf 
merger is defined to consist of the region primarily sup- 
ported by pressure, and not rotation (P < ^pw^), while 

the disk is rotationally-supported {^pv'^ > P). Further- 
more, our definition does not separate out the biconical 
jets and outflows into distinct additional components; 
the coronal region includes both the biconical jets and 
outflows. 
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2.2.2. Growth of the Magnetic Fteld and Magnetic Stress 

We expect the disk to be strongly unstable to the 
magnetorotational instability, a nd in deed, we confirm 
this to be the case. In figure 3(a) we show the de- 
velopment of the magnetic energy, for each run. For 
higher resolutions, we achieve a higher peak magnetic 
energy, though the runs exhibit a trend towards conver- 
gence in the peak magnetic energy with increased reso- 
lution. Specifically, the difference in the peak magnetic 
energy of run H-bh {E'p^^^ = 1.44 x 10*® erg) and S-bh 

{E'^i = 1-38 X 10*8 erg) is 0.46 that of the difference 
between the next two lowest-resolution models, S-bh and 
L-bh {E"^^^^ = 1.25 X lO-** erg). 

Next, we explore the role of stresses within our model. 
The r-0 component of the stress tensor. 



Tr4, = pSvuSv^ 



BrB^ 
47r 



(7) 



governs the transfer of angular momentum in the disk. 
The first term of the stress tensor is the Reynolds stress, 
and the second, the Maxwell stress. Here Svb. and Sv^ are 
the fluctuations of the radial and azimuthal velocity com- 
ponents. Analytically, during the linear growth phase of 
the MRI in a near-equilibrium disk, these fluctuations are 
the departures of the local fluid velocity from a circular 
orbit — specifically, Svr = vr and 5vci> = — i?r2(i?), 
where fl{R) is the disk angular velocity at radius R. 
However, it is well-known that due to the large turbulent 
fluctuations in fully-developed MRI turbulence, defining 
the mean disk velocity is problematic, even in a time- 
averaged sense, and consequently there is no unique pre- 
scription for specifying the Reynolds stress in a global 
disk s imulation — see for instance, |Hawley fc Krolik] 
(|2001|). Here, we focus upon the Maxwell stress as a 



proxy for the total stress, since 3D simulations of the 
MRI typically find the Maxwell stress dominates the 
Reynolds stress b y factors of 3- 6 — see, for instance, 
Davis et al. ( 2010 ). In figure 3(b)[ we plot the ratio of the 
spatial average of the Maxwell stress to the gas pressure, 
which defines an effective magnetic Shakura-Sunyaev pa- 
rameter am = {~ ^BuBfj,) I {Pgas), as a function of 
time, for each run. Here, brackets indicate spatially- 
averaged quantities over the disk and coronal regions, 
excluding the white dwarf merger itself. 

We find overall good agreement in the Maxwell stresses 
in all models computed, with all models apart from S-bm 
converging towards a value of ~ 0.01. Moreover, our 
standard model, which has been evolved for the longest 
runtime, has a relatively steady am ^ 0.01, indicative of 
sustained accretion and angular momentum transport. 

The white dwarf merger, which is initially weakly mag- 
netized with a mean, purely poloidal magnetic field of 
2.8 X 10*^ G, becomes rapidly magnetized over several in- 
ner rotational periods, as the MRI develops in the disk 
an d mag netic flux is advected into the merger — see fig- 



4(a) The field strength is time- variable, particularly 
at early times, which is expected given the turbulent non- 
steady nature of the MRI in the disk — just as the mass 
accretion is highly variable, so too must be the advec- 
tion of magnetic fiux into the merger. At late times, the 
total mean magnetic field strength within the merger is 
~ 2 X 10® G, typical of HMFWDs. The ratio of the 
mean toroidal field strength to the mean poloidal field 



strength is shown in figure [4(b)| The final field IS pre- 
dominantly toroidal, with a mean value of Bt/Bp ~ 1.5. 
Both the field strength and the geometry, as traced by 
the toroidal to poloidal field strength ratio, are highly 
time- variable, indicative of a disordered interior magnetic 
field. Given our limited resolution within the white dwarf 
merger of roughly 60 cells in our standard model S-bh, 
our final field strengths are very likely limited by numer- 
ical resistivity. Our results do, however, demonstrate 
that HMFWD field strengths may be achieved through 
the double-degenerate channel. 

The global magnetic energy ii^mag exhibits a slow decay 
at late times. About 7 x 10''^ erg, or roughly half of 
the total drop in magnetic energy is actually due to the 
outflow of magnetic fiux from the problem domain in 
our model S-bh. At late times, more magnetic fiux is 
outfluxed than the net magnetic energy generated within 
the domain by the MRI, as the MRI saturates and slowly 
winds down, and turbulent magnetic energy is dissipated 
through reconnection. The net result is a net decrease in 
magnetic field energy on the problem domain. 

2.2.3. Mass Accretton and Mass Outflow 

The stresses developed by the MRI drive mass accre- 
tion through the disk. Additionally, previous studies 
have demonstrated that the turbulence produced within 
the disk can develop coronal mass outflows (Machida] 
et al.|2000[|Miller fc Stone|2000[ [Suzuki fc Inutsuka,200' 



Flock et al. 201 1[ ). The implications of significant mass 
outflow are particularly profound for early-time optical. 
X-ray, and radio observations of SNe la, which are sensi- 
tive probes of the circumstellar environments surround- 
ing the progenitor systems. 

We characterize the initial mass of the white dwarf 
merger, disk, and magnetized corona by analyzing each 
of these at an early point in the simulation {t < 500 s), 
where the MRI is fully-developed, but prior to signifi- 
cant post-merger mass changes. At these early times, 
the white dwarf merger mass is 0.96 Mq by our specified 
criteria; the disk and the corona account for 0.20 Mq and 
0.04 Mq, respectively. 

In addition to the mass accretion from the disk onto 
the white dwarf merger, a significant fraction of the to- 
tal mass of the disk is lost, primarily through an outflow 
driven near the interface of the corona with the bicon- 
ical jet. Over the duration of our standard model, we 
determine that the disk has lost nearly 90% of its ini- 
tial mass over the duration of the simulation, through a 
combination of accretion and outflow. Of this amount, 
just over 82%, or 0.16 Mq, is accreted onto the white 
dwarf merger, with the remainder either transferred into 
the corona or outfluxed and lost from the domain. In 
total, nearly 0.06 Mq is outfluxed from the domain and 
subsequently lost to the simulation. 

2.2.4. Spin-Down of White Dwarf Merger 

A further key question is the possible influence of the 
magnetic field upon the spin of the white dwarf merger. 
Previous studies of core-collapse supernovae have sug- 
gested that the outward transport of angular momentum 
results in a spin-down of the p rotoneutron star and it s 
surroundings — see for instance, Thompson et al. ( 2005 ). 
In the context of the double-degenerate model of SJNe la, 
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Fig. 4. — a) The evolution of the mean magnetic field within the white dwarf merger versus time, shown for both the mean toroidal 
(solid) and mean poloidal (dashed) lines for model S-bh. b) The ratio of mean toroidal to mean poloidal magnetic field within the merger, 
versus time, again for model S-bh. 



a loss of rotational support of the central white dwarf 
merger may yield additional or perhaps even dominant 
compression beyond that provided by accretion from the 
disk. 

We find the central white dwarf merger is spun down 
significantly on a very rapid timescale, due primarily 
to the development of Maxwell stresses at the bound- 
ary of the white dwarf merger. In particular, we find 
that the merger, whose initial angular momentum is 
2.6 X 10^° g cm^ s""'^, spins down significantly to 8.0 x 
2^q49 g (,j-|^3 g-i ^Yie end of the standard model simu- 
lation. This amounts to a brak ing timescale of J/\J\ — 
several xlO^ s. In figure 5 (a 

averaged angular velocity il as a function of cylindrical 
radius r. The white dwarf merger is seen to be in solid- 
body rotation at r < 10® cm, with the outer portion of 
the domain is in Keplerian rotation. We note the taper- 
ing of the angular velocity plot at radii r > 6 x 10® cm 
at t = s is an artifact of the initial SPH particle distri- 
bution, which lacked any particles in this region. Addi- 
tionally we note that the shear in the innermost portion 
of the disk, near r = 2 x 10® cm, flattens out as the free 
energy available in the shear is tapped to drive the MRI. 
As time evolves, the merger is seen to spin-down signif- 
icantly, from n ~ 0.3 initially to almost half of that 
value, ^ 0.17 s^^, finally, at the end of the run. 

In order to verify that the simulated spin-down of the 
merger is indeed physically accurate, and not the conse- 
quence of numerical errors, we have computed the spin- 
down expected from the conservation of angular momen- 
tum and compared this to the measured spin-down. In 
particular, in axisymmetry, the inviscid angular momen- 
tum evolution equation can be written as 



represents the divergence of the angular momentum flux. 



R 



(9) 



we plot the vertically- ^^^^^ gg^^^ ^ 



In post-processing, we integrate the divergence of the 
angular momentum flux F^f, inside the merger, and com- 
pare it with the actual angular momentum evolution of 
the merger as obtained by the full simulation. For the 
purposes of this computation we define the merger as the 
fixed region with a distance less than 1.5 x 10'' km from 
the domain center, which is consistent with our previ- 
ous definition base d upo n pressure support. The result 



d_ 
di 



(pRv^) + V-R 



= 



(8) 



Here, Bp is the poloidal magnetic field. The second term 



The cumulative effect of the 
magnetic and hydrodynamic stresses are integrated over 
volume and time to produce a net change in angular mo- 
mentum, shown in the dash-dotted curve. The Maxwell 
stresses are roughly three times more significant than the 
Reyno lds stress, simila r to previoiis MRI studies, e.g.. 
Flock et al. (20ll| and [Davis et^ ( |2010[ ). 

'I'he sum of these predicted net changes in angular 
momenta are plotted, and added to the initial angular 
momentum, as the dashed curve, which compares with 
the actual simulated angular momentum, as shown in 
the solid curve. The results agree to within 10%, which 
establishes that the spin-down is physically driven by 
MRI-generated Maxwell stresses, and is not the result 
of numerical or artificial viscosity. Furthermore, by sub- 
sampling the output interval, we have determined the 
dominant error in the calculation to be the finite output 
cadence of 10 s. Thus, while our calculation places the 
upper-bound of the numerical errors at ^ 10%, the true 
numerical error is very likely to be much smaller. 

2.2.5. Nuclear Burning 

The question of whether the merger will initiate a det- 
onation, and possibly result in a SNe la, hinges crucially 
upon the peak nuclear burning rate. Because the equal- 
mass merger case produces a central peak in tempera- 
ture, the burning rate achieves its highest value at the 
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(a) Q vs. r (b) J, AJ vs. time 

Fig. 5. — a) A plot of the angular velocity Q versus cylindrical radius r, shown at various times — solid curve (t = s), dashed curve 
[t = 10* s), and dotted curve {t = 2 X 10* s) for model S-bh. b) A plot of the angular momentum and computed changes in angular 
momentum of the white dwarf merger vs. time for model S-bh. The solid line represents the actual simulated spin, whereas the dashed 
line is the computed spin based upon the combined effects of the Maxwell and Reynolds stresses, whose cumulative effect is shown in the 
dash-dotted and dotted curves, respectively. 



merger center. In figure 2.2.5 we plot the central con- 
ditions for the white dwarf merger in the log T — log p 
plane, as a thick solid line. Also shown are equicon- 
tours of both the carbon-burning (dot-dashed) and neu- 
trino timescales (dashed), at values of 10^ yr, 10^ yr, and 
10~^ yr. The critical condition for thermal runaway, at 
which the timescale for carbon burning is equal to that 
of neutrino cooling, is shown as a thin solid line. The de- 
generacy boundary in the temperature-density plane is 
demarcated by the dotted where the temperature equals 
the Fermi temperature: T = Tp. The central condi- 
tions of the merger remain partially, though not fully, 
degenerate throughout the evolution. At the beginning 
of the simulation, at t = s immediately after the ini- 
tial merger, the neutrino cooling time is shorter than 
the carbon-burning time at the center of the merger. 
As the white dwarf spins down and accretes mass from 
the disk, the central region compresses further, and the 
central conditions become thermally unstable, with the 
timescale for carbon burning shorter than neutrino cool- 
ing at around ^cc = 10^ yr. At this point, the central 
region of the white dwarf has ignited. 

After ignition, the central temperature of the merger 
continues to rise, reaching a peak temperature T ~ 
10^ K. The burning timescale at this temperature is 
roughly 1 yr, which is much longer that what we can 
feasibly follow in this set of simulations. We return to 
this issue below in the conclusions. 

3. COMPARISON WITH PREVIOUS RESULTS 

While there is a large body of work o n the final stages 
of the merger of white dwarf bina ries (Mochkov itch fc 
Liviolll989( Tl990l IRasio fc Shapiro|fT 995, .Segretain et a 



1997; Guerrero ct al. 2004; 'D^ouz a'StaT 
et al.,,2007. .Motl et al..,2007. . Pakmor et aLn^w^^, 
et al.|[20TT Zhu et al. 201 1[ ) there are relatively few in- 



|2006| lYoon 
| |2010( |Dan 



depth studies of the post-merger phase of evolutio i i. Sey - 



Saio fc Nomotoj ( |1998 ), and |Saio fc Nomoto| ( |2004[ ), 
follow the evolution oi a spherical remnant accreting 
near the Eddington limit, and concluded that the off- 
centered ignitio n would lead t o an a ccretion-induced col- 
lapse. Later, Yoon et al. (20071 followed six white 



dwarf mergers using SFH, varying both the mass ra- 
tio and the number of particles. Their most-advanced 
model consisted of a super-Chandrasekhar SPH simula- 
tion of 0.6 Mq -|- 0.9 M0 white dwarf binary, which was 
followed through merger and a short time (300 s) af- 
ter merger using invisc id hydrodynam i cs. A dditionally, 



Shen et al.| ( |2012l ) and Schwab et aL] ( |2012[ ) follow the 
post-merger evolution 01 a range of super-Chandrasekhar 
mass white dwarf binaries in ID spherical and 2D spher- 
ical geometry, respectively, using a Shakura-Sunyaev vis- 
cosity prescription for the transport of angular momen- 
tum. The comparison of our work to these previous mod- 
els is complicated by our focus upon an equal-mass sub- 
Chandrasekhar case, whereas most prior models in these 
studies focus upon unequal mass, super-Chandrasekhar 
models. We are, however, able to compare the broadest 
features and underlying physics of all models. 
The morp holog y of the merge r s stud ied by both |Yoon| 



et al. (|2007|) and Schwab et al. 

wart core 



(20121, which consist of 



lot envelope, and an ac- 



a rotating white dwart core, a 
cretion disk is in broad agreement with the structure 
of our initial SPH conditions, with the crucial difference 
that their unequal mass models have temperatures which 
peak off-center, as other stu dies have also shown — e.g. 
Loren-Aguilar et al.l ( |2009l ) and |Zhu et al.| ( [20TT| ). We 



do note that our initial central peak teinperatu r es are 



somewhat higher than those found by Zhu et al 
for similar sub-Chandrasekhar mergers, and that 



poll I 

this is 



likel y due to d ifferences in the SPH modeling. Addition- 
ally, Zhu et al. ( [2011) find that peak central temperatures 
are achieved only for non-synchronously rotating white 
dwarfs, as ours are. 



eral classic spherical studies of Saio fc Nomoto (1985), 



Yoon et al. (2007) parameterize the end state of their 
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Fig. 6. — A phase plane of the central conditions of the simulated white dwarf merger for model S-bh. See text for description. 



SPH models by mass, and advance a range of models 
in a ID hydrodynamic code incorporating a description 
of angular momentum transport via hydrodynamic in- 
stabilities with an imposed timescale of 10^ — 10^ yr. 
Our models, which include self-consistent angular mo- 
mentum transport through th e magnetorotational insta- 
bility, as well as the models of Schwab et al. (2012) using 
an a model, find that the tiniescaie tor accretion of the 
disk is consistent with the turbulent viscous accretion 
timescale, and many orders of magnitud e shorter than 



model assumptions of Yoon et al. 
are, however, consistent wit h the 



(20071. Our findings 



imescale of turbulent 



disk accretion suggested by van Kerkwijk et al. (2010). 
Our findings are also consistent with recent experiments 
which find that hydrodynamic tur bulence is inefficient at 
transporting angular momentum (Ji et al. 2006). They 
are gener ally incon sistent with accret ion at the Edding- 

ton rate ( |Saio fc No moto 1985, 1998', |2004| . 

Furthermore, SEen et al., (,2012) and [Schwab et aL 



(2012) find that the influence of a spatially-dependent 
turbulent viscosity (a ~ 0.03) drives their mergers to 
spin down completely over a timescale of ~ 3 x 10^ s, 
quite similar to our finding with the full MRI that the 
white dwarf merger is magnetically braked on a compa- 
rable timescale. Moreover, they further find that as tur- 
bulent viscosity dissipates rotational energy into heat, a 
thermally-supported outer region develops , also broadly 
similar to our heated corona. Additionally, [Schwab et al.| 
(2012) find evidence for sustained carbon burning, as we 
do in our models. 

However, there is a subtle but important distinction 
in the physics of the angular momentum transport and 



heating mechanisms at work here. The Shakura-Sunyaev 
prescription posits that shear energy is locally dissipated 
as heat energy over a viscous timescale. In contrast, in a 
magnetized system, the conversion of shear energy in the 
disk does not directly lead to heating. Instead, under the 
action of the MRI, shear energy is converted into mag- 
netic energy, which is subsequently buoyantly displaced 
into the corona. The corona is ultimately heated through 
the non-local dissipation of magnetic energy into heat en- 
ergy. In our models, this is accomplished through numer- 
ical resistivity, though in reality the magnetic dissipation 
will occur through both physical resistivity and rapid re- 
connection. Similarly, the inner white dwarf merger is 
heated under compressional work from both accretion 
and spin-down. The end result in our case is a highly 
inter mitten t and non-uniform thermal structure — see 
Indeed, about 7 x 10'''' erg, or roughly half 
magnetic energy, is out fluxed entirely from 



figure 1(b)] 
of our peak 

our problem domain. In contrast, [Schwab et al. (2012) 
find their outer envelope is smoothly heated through the 
action of turbulent viscosity, which is the direct result of 
the conversion of local shear energy into heat. Further 
magnetized models will be able to determine the precise 
extent to which the thermal structure of the hydrody- 
namic and magnetic models differ, and identify possible 
implications for nuclear burning and detonation. 
The most s ignific ant difference, however, between 



Schwab et al. (2012) an d the models p resented here re- 
late to mass oiattluxes. Schwab et al. ( [2012^ ) find that 
mass outflows are limited to less than 1O~^M0 in their 
standard 0.6 Af0 -I- 0.9 model. In contrast, our sub- 
Chandrasekhar model drives a vigorous outflow several 
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orders of magnitude greater. F| Our results are consis- 
tent with global MRI studie s, which often find signif icant 
outflows — s ee for instance, De Villiers et al^ ( 2005 ) and 



Flock et al. (2011). ine outflows generatea m magne- 
tized disics are fundamentally driven by the buoyancy 
of the magnetized corona and its interaction with the 
strongly-magnetized bi conical jet. The Shak ura-Sunyaev 
prescription utilized by Schwab et al. (2012) captures the 
basic aspects of mass and angular momentum transport 
by specifying the r-(f) and 9-(f) components of the stress 
tensor, but lacks both magnetic pressure and tension and 
the Lorentz force, which play a crucial role in driving 
magnetized outflows. 

Moreover, while our models suggest that 
magnetorotationally-driven outfluxes may influence 
the immediate circumstellar environment surrounding 
the white dwarf merger, the precise fate of this outfluxcd 
matter is underdetermined in our simulations, which do 
not capture the sonic surface {Rs ~ 2GM/cl ~ 10^^ cm 
for T ^ 5 X 10^ K gas) of the outflowing material 
within the simulation domain. There is some evidence 
based upon visualizations of the density and magnetic 
field that our result may be more complex than a 
pure infall or thermal wind outflow, and may consist 
of a combinat ion of a thermal wind and reconnection- 
^ The 

net 



driven infall ( Igumenshchev et al 
outflow 



course. 



total 
central to 



mass outflow and mtall is, of 
the determination of the nuclear energetic yield and 
corresponding light curves of any possible SNe la that 
originates from a double-degenerate merger — not 
only of sub-Chandrasekhar mass as studied here, but 
potentially also of super-Chandrasekhar mass if the 
merger does not result in a prompt detonation. 

Furthermore, the topology of the seed field plays a 
cruci al role in determining whether a large-scale ordered 
field (Beckwith et al. 2008). Consequently, our findings 
with regard to magnetic braking, as well as the jet and 
outflow, are to some extent dependent upon our initial 
seed magnetic field, which consists of a single torus of 
magnetic flux. Further work will need to be done to de- 
termine the degree to which these conclusions are robust 
in the presence of realistic seed magnetic fields generated 
from the white dwarf merger process itself. 

4. CONCLUSIONS 

The nuclear burning time at the endpoint of our stan- 
dard simulation ^cc 1 yr is much longer than our 
evolutionary time of 6 hours. Thus, while the central 
temperature continues to rise even throughout the sim- 
ulation, the outcome of the nuclear burning is not yet 
clear. The conditions are, however, favorable to super- 
sonic detonations, as opposed to subsonic deflagrations. 
At a density ^ 10^ g cm"-^, the specific energy release 
in carbon burning exceeds the specific internal energy by 
about an order of magnitude, which gives rise in an over- 
pressure also about an order of magnitude gre ater than 
the ba ckground pressure — see, for instance, |Nomoto| 
(1982). Thus our results suggest that under appropriate 



^However, in an unpublished synchronously-rotating, equal- 
mass, super-Chandrasekhar 0.75 Af0+O.75 Mq run, they do obtain 
a thermally-driven expansion of their white dwarf merger, which 
would amount to 0.1 A/0 beyond our domain (J. Schwab, private 
communication). This is similar to the total amount of outflux we 
see in our lower-mass, sub-Chandrasekhar model. 



conditions of near-equal mass non-synchronous mergers, 
a sub-Chakdrasekhar merger may give rise to a central 
detonation powering a SNe la. However, while this sce- 
nario is promising, a detonation is not assured, as the 
detonation initiatio n condition depe nds on the tempera- 



ture profile (Seitenzahl et al._ 
ther computational studies wi. 



20091). Consequently, fur- 
1 need to be carried out to 



explore the outcome in greater detail. 

Additionally, while the 2D axisymmctric models pre- 
sented here begin to shed light on the role of the mag- 
netic field in binary white dwarf mergers, both the accre- 
tion itself and the development of the Maxwell stresses 
leading to magnetic braking of the rotating white dwarf 
are necessarily limited in 2D axisymmetry. Additionally, 
some of the long-lived magnetic flux tori which we see 
in 2D may or may not prove to be stable in 3D, which 
may significantly alter the evolution of the corona. A 
fuller picture of the evolution of the rotating white dwarf 
merger will require 3D simulations, and while these will 
be demanding, we expect that these will should be more 
favorable to thermonuclear runaway. 

An off-centered ignition, possibly leading to an 
accretion-induced collapse, may still be possible for sig- 
nificantly unequal mass mergers, but we expect that the 
infiuence of the magnetic field will alter some conclu- 
sions of previous models quantitatively, if not qualita- 
tively. Furthermore, we expect that for unequal mass or 
significantly lower-mass mergers, the peak temperatures 
reached even during the post-merger phase may never 
reach ignition conditions. In these cases, the outcome 
will not be either a SNe la or an accretion-induced col- 
lapse, but rather a stable object. Our findings are con- 
sistent with theoretical models which have predicted the 
growth o f the white dwarf magnetic field throug h binary 



systems (|Tout et al 



Nordhaus et al. 2011 



Garcfa- 



Berro et "aT ' 2U12p . Our simulations demonstrate that 
such an object will be strongly-magnetized, and may ac- 
count for the existence of HFMWDs. Moreover, if both 
the double-degen erate channel of wh i te dw arf mergers, 
as conjectured by Garcia-Berro et al. (2012), and a sep- 
arate channel o f white dwarf-low mass companions, as 
conjectured by Nordhaus et al. ( 2011|) , are simultane- 
ously active, we expect the Hl^'MWD distribution to be 
bimodal in field strength. 

HFMWDs and SNe la may therefore represent dis- 
parate outcomes of the same fundamental astrophysi- 
cal process of merging white dwarfs. Consequently, the 
HMWWD magnetic field distribution may help inform 
our understanding of which double-degenerate systems 
are failed SNe la, yielding stable white dwarfs as op- 
posed to thermonuclear detonations. Further magnetized 
models of a wide range of white dwarf masses, rang- 
ing from ONe white dwarfs to He white dwarfs through 
low-mass stellar and planetary companions, will help es- 
tablish the HMFWD birth field distribution, which can 
then be modeled to comp are directly against their ob - 
served values in the field ( Vanlandingham et al. 2005). 
Such a study will also help more fully elucidate the condi- 
tions for thermonuclear runaway to SNe la in the double- 
degenerate channel. 

As we have demonstrated, the magnetorotationally- 
driven disk turbulence produces a significant outflux of 
mass from the disk. Although this effect has not been 
pointed out in the context of merging white dwarfs and 
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SNe la to the best of our knowledge, we expect such out- 
fluxes to be a general outcome of the MRI during the 
post-merger stage of evolution of the double-degenerate 
system. Therefore, our models suggest that the immedi- 
ate circumstellar environment surrounding white dwarf 
mergers may not be as clean as previously believed. To 
date, only upper limits for the X-ray luminosity for both 
53 SNe la using Swift ( [Russell fc lm mler"201 2l), as well 
in the X-ray and rad io tor 2Ullte (fHorcsh et al. 

have been established 



[Chomiuk et al.||2012 ) 
magnetized models or 



2012 
Further 



magnetized models on extended domain sizes. 



using a 

range of seed fields, and capturing both the explosion 



and the sonic point of the outflow material may help 
elucidate a precise, predictive observational signature of 
the double-degenerate channel of SNe la. Such simu- 
lations will address the final fate of the outflows — a 
wind, infall, or a combination thereof. They will also 
address to what extent if any the outfluxes may have on 
the light curve and spectra at early times, as the super- 
nova blast wave encounters outfluxed circumstellar mat- 
ter surrounding the double-degenerate system, or infall 
back onto the white dwarf merger powers the light curve 
through accretion energy. 



APPENDIX 



In this appendix, we detail the azimuthal-averaging procedure used to average the 3D SPH particles onto an axisym- 
metric Eulerian r-z mesh. A simple method involves averaging the 3D SPH particles onto a 3D cylindrical Eulerian 
mesh, then angle-averaging this mesh. However, a more elegant and efficient method, which we adopt, is to break 
up the angle-dependent quantities, which involve only the smoothing kernel, from the particle data. This reduces the 
angle-averaging procedure to simple quadrature, with the weights written as integrals to be pre-tabulated once. The 
averaged quantities are then rapidly and efficiently calculated by weighting the particle data at each point on the r-z 
mesh over lookups of the pre-tabulated weights. 

The kernel function is the central mathematical object within the cont ext of SPH which allows cont inuous fluid 
quantities to be calculated from discrete particle data. We adopt the form ( Monaghan fc Lattanzio|[l985 ) : 



W{Y,h) 




2(|-i/2), 0<z.<i 



Here v is the dimensionless distance scaled to the smoothing length 
dimensionless kernel function W as 



(1) 



\v\/h. We further define a rescaled 



W{v) = -7rh^ W{r,h) 



8 (3 
3 U 

3 ^2 

0, 



(2) 



W depends solely on i/. In the "scatter" interpretation of SPH, each quantity A at a specific position r is given by 



N 
i=l 



W{\r-ri\,h,), 



(3) 



with the number of particles N and subscripts denoting quantities for a specific particle. To angle-average this quantity 
A, we integrate it over the angle if and divide it by 2 tt: 



A{r, z) 



1 

2^ 



2Tr 



dLpA{r,(p,z) 



(4) 



Because the quantities m^, pi and Ai do not depend upon angle, we may interchange the order of integration and 
summation, and write A{r^ z) as 



A{r, z) 



1 



N 



1 



(5) 



Written in this way, we can transparently see that in the scatter interpretation, angle-averaging involves only the 
smoothing kernel itself. We can further simplify this expression by defining the dimensionless 2D cylindrical distance 
— ■\/(?' — r^)^ -\- {z — Zi)^/hi and dimensionless ratio x — 2rri/hf to express 



u = V '^2D^ + x[l - cos{ip - ipi)] (6) 
Next, we define Wo{i'2Dt x) to be the angle-averaged dimensionless kernel, expressed as a function of V2D and x: 

3 



Wo(i^2D,a;) 



87r2 



dipW{l^2D,X,C0s{ip)) 



(7) 
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This allows us to efficiently compute the azimuthally-averaged temperature T and the z-velocity Vz at each point as 



N 



P{r,z) =^^Wo{v2Yi,x), T{r,z) = ^Ti^^Wo{v2T>,x), and Vz{r,z)=^v, 



N 



N 



i=l 1=1 i=l 

One can similarly show that the r and (j) velocity components can be expressed as 

N 

Vrir,z) = 



and 



where we have defined 



and 



i=l 



AT 



TTi ' TTl " 



Wo{iy2B,x). (8) 



(9) 



,sinv^2D; 



Wo 



Wo ,sin(,^2Di 



dipW{v2T>TX,ip) COS ip 



87r2 



27r 



d(/3 VF(z^2D, a;, v) saiLp 



(10) 



(11) 



(12) 



analogously to Wq(i>2T),x). In practice, the functions Wo(j^2D,a;), Wq, 005(^^20, a;), and Wo,sin(j^2D, a;) can be pre- 
tabulated as functions of both V2D and x. In this way, one can rapidly and efficiently calculate Eulerian grid quantities 
from azimuthally-averaged SPH data. 



Simulations at UMass Dartmouth were performed on a computer cluster supported by NSF grant CNS-0959382 and 
AFOSR DURIP grant FA9550-10- 1-0354. PT, GJ, and DL acknowledge support at the University of Chicago in part by 
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